Low-energy spin dynamics in the [YPc 2 ]° 5 = 1/2 antiferromagnetic chain 
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1 H nuclear magnetic resonance (NMR) measurements in [YPC2] , an organic compound formed 
by radicals stacking along chains, are presented. The temperature dependence of the macroscopic 
susceptibility, of the NMR shift and of the spin-lattice relaxation rate 1/Ti indicate that the unpaired 
electron spins are not delocalized but rather form a S = 1/2 antiferromagnetic chain. The exchange 
couplings estimated from those measurements are all in quantitative agreement. The low-energy 
spin dynamics can be described in terms of diffusive processes and the temperature dependence of 
the corresponding diffusion constant suggests that a spin-gap around 1 K might be present in this 
compound. 

PACS numbers: 76.60.Es, 75.10. Pq, 75.40.Gb, 75.10.Jm 



I. INTRODUCTION 

Molecular solids have attracted much interest since 
decades owing to the possibility to easily tune their prop- 
erties either through a chemical bottom up approach or 
by varying physical parameters as the external pressured 
One of the most versatile families of molecular solids 
is the one based on phthalocyaninc (Pc= C^H^Ns) 
molecules.— In fact, the employment of these materials 
in different areas, ranging from the fabrication of organic 
light emitting diodes, to contrast agents or spintronics 
materials, has been envisaged. Pc-based compounds have 
attracted significant interest in the last decade after it has 
been suggested that high temperature superconductivity 
could be induced in these materials by alkali doping 3 and, 
more recently, when it has been recognized that neutral 
[LnPc2]° molecules, with Ln a lanthanide ion, are molec- 
ular nanomagnets with extremely long coherence times at 
liquid nitrogen temperature^— Owing to the flat shape of 
Pc molecules, the structure of Pc-based materials is typ- 
ically characterized by chains along which Pc molecules 
tend to stacks Accordingly some of the Pc-based mate- 
rials show many similarities to the Beechgaard salts J£ 

Bis(phthalocyaninato) yttrium [YPC2] compound can 
be considered the parent compound of the aforemen- 
tioned [LnPc2]° molecular magnets. In fact, it is charac- 
terized by the absence of localized / electrons and the 
microscopic properties are mainly associated with the 
presence of an unpaired electron delocalized in the 02 tt 
orbital, due to the one-electron oxidation of the [YPc 2 ] ~ 
unitiii Thus, [YPC2] allows to investigate the spin dy- 
namics associated only with this unpaired electron spin, 
independently from the one due to / electrons. One of 
the most suitable tools to address this aspect is nuclear 
magnetic resonance (NMR) technique. In this work we 
present an experimental study of the magnetic proper- 
ties of [YPC2] compound by means of magnetization and 
nuclear magnetic resonance (NMR) measurements. The 
temperature dependence of the macroscopic susceptibil- 
ity, of the NMR shift and of the spin-lattice relaxation 
rate 1/Ti clearly show that this system is a prototype 



of a S — 1/2 antiferromagnetic chain, characterized by a 
diffusive low-frequency spin dynamics and, possibly, by 
the presence of a low-energy spin-gap. 
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FIG. 1: Temperature dependence of static uniform suscepti- 
bility xs for [YPC2] complex, derived from SQUID magneti- 
zation measurements. The solid line shows the best fit of the 
data to Curie- Weiss law. 



II. EXPERIMENTAL RESULTS AND 
DISCUSSION 

[YPC2] polycristalline samples were synthesized by us- 
ing some modifications of the protocol published in Ref. 
|12| . All reagents were purchased from Across or Aldrich 
and used without further purification. A mixture of 1,2- 
dicyanobenzene (42 mmol), Y(acac)3 -4H20 (5 mmol), 
and l,8-diazabicyclo[5,4,0]undec-7-ene (DBU) (21 mmol) 
in 50 mL of 1-pcntanol was refluxed for 1.5 days. The 
solution was allowed to cool to room temperature. The 
precipitate was collected by filtration and washed with 
rt-hexane and Et20. The crude purple product was re- 
dissolved in 800 ml of CHCl 3 /MeOH (1/1) and undis- 
solved PCH2 was filtered off. Both forms, blue (anionic 
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[YPc2] - ) and green (neutral [YPC2] ), were obtained si- 
multaneously, as revealed by electronic absorption spec- 
tra. In order to convert the unstabilized anionic form to 
the neutral one, the reaction mixture was presorbed on 
active (H2O-0%) basic alumina oxide. Purification was 
carried out by column chromatography on basic alumina 
oxide (deactivated with 4.6% H2O, level IV) with chloro- 
form methanol mixture (10:1) as eluent. In general, the 
yield was 30-35%. According to microelemental analysis 
based on atomic spectroscopic methods (ICP) performed 
at Mikroanalytisches Labor Pascher, the powder sample 
contains molecules of [YPC2] , water and CH2CI2 in ra- 
tio 1:1:1/3. The molecules crystallized in the space group 
P2i2i2i (7-phase), as reported in Ref. O 

DC magnetization (M) measurements have been per- 
formed by using an MPMS-XL7 Quantum Design super- 
conducting quantum interference device (SQUID) mag- 
netometer. The magnetization was found to depend 
linearly on the magnetic field intensity H, for H < 5 
kGauss, over all the explored temperature range and, ac- 
cordingly, the macroscopic static uniform susceptibility 
can be written as xs = M / H. The temperature de- 
pendence of xs reveals the presence of antiferromagnetic 
correlations. In fact, xs(T) can be nicely reproduced by 
a Curie- Weiss (CW) law 
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where C = g 2 ii 2 B S{S + 1)Na/ {Sks) is Curie constant (/xg 
the Bohr magneton, g the Lande factor, Na Avogadro's 
number and ks Boltzmann constant). is the CW tem- 
perature and xo a temperature independent term mainly 
due to diamagnetic and Van-Vleck corrections. By fitting 
the data, leaving all three parameters free, we found an 
antiferromagnetic CW temperature = 5.37 ± 0.04 K 
(Fig. [J) and a Curie constant C = 0.342 ± 0.002 erg 
• K/G 2 , quite close to the value 0.375 erg-K/G 2 , ex- 
pected for a S = 1/2 system. If we fixed C = 0.375 
erg-K/G 2 the fit was still good and the CW temperature 

= 6.18 ± 0.03 K. The temperature dependence of xs 
shows that the unpaired electron spins are localized along 
the chains formed by [YPC2] molecules, although a cer- 
tain overlap of the it orbitals of adjacent molecules must 
be present in order to justify the magnitude of the antifer- 
romagnetic exchange coupling J e . In fact, although this 
system should present a narrow half-filled band, the size- 
able Hubbard Coulomb repulsion U ~ leV prevents the 
electron derealization along the chainJ^ In this limit, 
J e = = 4t 2 /U, with t the hopping integral among 
adjacent molecules. From the estimated value of one 
would derive a band width formed by the overlap of 02 
orbitals in adjacent molecules W = At ~ 0.05 eV<C U, 
justifying the spin localization along the chain. 

The X H NMR spectra were obtained in the 1.6-300 K 
temperature range for magnetic field intensities H = 9 T, 

1 T and 0.3 T. The spectra were derived from the Fourier 
transform of half of the echo formed after a7r/2 — t — tt 
pulse sequence, when the full NMR line could be irradi- 
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FIG. 2: Temperature dependence of 1 H full NMR linewidth at 
half intensity in [YPC2] , at three different magnetic fields. In 
the inset the linewidth is normalized by the Larmor frequency 
after subtracting the constant linewidth due to nuclear dipole- 
dipole interaction. 



ated or, otherwise, from the envelope of the echo ampli- 
tude upon varying the irradiation frequency. The line- 
shape was gaussian in all the investigated temperature 
range. For H — 9 T and H = 1 T a broadening of the 
spectrum can be observed at low temperature, which is 
more pronounced at higher field intensities (Fig. [5]). On 
the other hand, for H = 0.3 T the linewidth Av is prac- 
tically temperature independent and the broadening is 
likely to be due just to nuclear dipolc-dipole interaction. 
The increase of the linewidth with H suggests that the 
low temperature line broadening originates from some 
anisotropy in the hyperfine coupling, which for a powder 
gives rise to a linewidth proportional to the susceptibility. 
In fact, it is noticed that if we subtract the T-independent 
contribution at H = 0.3 Tesla from the raw data and di- 
vide the linewidth by the Larmor frequency Vq, the data 
at different fields overlap (inset to FigJ^J). This result also 
indicates that there is not an additional internal field due 
to the onset of a long-range magnetic order down to 1.6 
K. 

The NMR paramagnetic shift AK = {vr — v^/vo, with 
vr the resonance frequency, shows a more pronounced 
temperature dependence (Fig. [3]). As expected, it was 
found to increase upon cooling, according to 
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namely the temperature dependence of AK should be 
the same of the macroscopic spin susceptibility. In fact, 
also AK(T) is found to obey a Curie- Weiss law with a 
Curie- Weiss temperature = 7 AK ± 0.3 K, close to the 
one derived from SQUID magnetization measurements. 
The small difference between those two type of measure- 



3 



0.0004 - 




1 0.0003 - 

m U 

T < 

\ 0.0002 - 
L 0.0001 - 
\ 0.0000 - 


/ 






00 0.005 0.010 0.01 5 0.020 0.025 

X s (emu/mole) 


|H - 9 T | 





50 100 150 200 250 



T(K) 

FIG. 3: Temperature dependence of 1 H paramagnetic shift 
AK in [YPC2] . The solid line shows the best fit of the data 
to Curie- Weiss law. In the inset AK is reported as a function 
of the macroscopic susceptibility. The solid line is the best fit 
according for a linear dependence of AK vs xs- 

merits could be associated with a tiny amount of im- 
purities which might affect the macroscopic susceptibil- 
ity. Accordingly, the measurement of the microscopic 
susceptibility with paramagnetic shift measurements is 
expected to provide a more reliable estimate of the static 
uniform spin susceptibility and of the Curie- Weiss tem- 
perature. By plotting AK as a function of \S a linear 
trend is attained (Fig. [3]) and from the slope it is possible 
to estimate the isotropic term of the hyperfine coupling 
tensor A = 180 ± 10 Gauss. 

The -"-H nuclear spin-lattice relaxation rate 1/Ti was 
measured in the 1.6 - 300 K temperature range and for 
different values of the external field. 1/Ti was extracted 
from the recovery of nuclear magnetization after a sat- 
uration recovery pulse sequence. The recovery law was 
found to be a single exponential in all the explored tem- 
perature range (Fig. 31 inset). This result is an evidence 
that the unpaired electron is delocalized onto a ir orbital 
within the molecule. In fact, since in the two phthalo- 
cyanine rings a large number of inequivalent proton sites 
is present, if the electron was on a more localized orbital 
a distribution of hyperfine couplings would be present 
and, accordingly, a stretched exponential recovery law 
should be observed. Moreover, the fact that hyperfine 
coupling seems quite isotropic indicates that it could orig- 
inate from the contact interaction between the unpaired 
electron spin in the 02 7r orbital and the H nuclei. 

The temperature dependence of 1 /T\ at different mag- 
netic fields is shown in Fig. (QJ. In general, for a relax- 
ation process driven by electron spin fluctuations one can 
write 

^ = ^£(IAJX«(q*<>U , (3) 
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FIG. 4: H nuclear spin-lattice relaxation rate temperature 
dependence for [YPC2] compound measured for different val- 
ues of the applied field. In the inset the recovery of the nu- 
clear magnetization as a function of the delay r between the 
saturating and the echo readout sequences is shown at two 
different temperatures. The solid lines show the best fit for a 
single exponential recovery. 

where 7 is the nuclear gyromagnetic ratio, \Aq\ 2 the form 
factor describing the hyperfine coupling with spin excita- 
tions at wave-vector q and S a ^ a (q, ujl) (a = x,y,z) the 
component of the dynamical structure factor at the Lar- 
mor frequency. In the high temperature limit, namely 
when the thermal energy is much larger than the ex- 
change energy (T > 6), the 1/Ti of a spin 5" = 1/2 
antiferromagnet becomes temperature and field indepen- 
dent and is given by-i£ 

+ , (4) 

where A x ~ A y ~ A are the components of the hyper- 
fine coupling tensor which is basically isotropic, while 
luh = (J e kB /h)y / 2zS(S + l)/3 is the Hcisenberg ex- 
change frequency, with z — 2 the number of nearest 
neighbour spins along the chain. By taking the mea- 
sured value of 1/Ti — 20 s _1 at high temperature, from 
Eq. dU it is possible to estimate the exchange frequency 
ll>h — 9.2 • 10 11 rad/s, corresponding to an exchange cou- 
pling constant J e ~ 7.0 K, in quite good agreement with 
the value which can be estimated from the NMR shift 
measurements. Upon decreasing the temperature, for 200 
K> T > 30 K, one observes a progressive slow increase 
of 1/Ti (Fig. (gj). In particular, it is noticed that nu- 
clear spin-lattice relaxation rate increases on decreasing 
temperature according to 

l/T x <xln x ' 2 {T /T) . (5) 

In fact, in Fig. (JSJ) one observes that (1/Ti) 2 is a linear 
function of 1/T, when reported in logarithmic scale. Re- 
markably, this logarithmic increase of 1 jT\ is expected in 
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a S = 1/2 Heisenberg antiferromagnet, but for T < Jj£. 
Here it is not clear why the logarithmic behavior extends 
up to T ^> J e . 
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FIG. 5: The 1/Ti squared is plotted as a function of T _1 , in 
logarithmic scale, for two values of the external field (H = 9 
T, circles and H = 1 T, squares). The dashed lines represent 
the best fits to Eq. ©. 

At about 20 K a peak in the nuclear spin-lattice re- 
laxation rate appears (Fig. whose intensity decreases 
by increasing the external field intensity. Eventually, be- 
low T ~ 5 K, the 1/Ti is only weakly temperature de- 
pendent. The maximum in 1/Ti, not associated with 
molecular motions, could be due to a form factor, which 
partially filters out the antiferromagnetic fluctuations as 
the system gets more and more correlated. 
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FIG. 6: The X H spin-lattice relaxation time 1/Ti in [YPc 2 ]° 
is plotted as a function of (oJo/2n)~ 1 ^ 2 for different selected 
temperatures. The solid lines show the best fit according to 
Eqs[6]and[7|in the text. 

The magnetic field dependence of 1/Ti (Fig. [5]) can 



originate from the diffusive nature of the spin correlation 
function, which in one dimension is characterized by long- 
time tails yielding to a divergence of the low-frequency 
spectral density J(ui)JZ In fact, in the presence of diffu- 
sive processes for the spin excitations 1 /T± can be written 
in terms of the spectral density for the spin excitations 
according to the following equatior 



,18. 



_J_ _ T k B Txo 
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\A 2 d J{oj ) +(A 2 + 7 -AYj J(oj e ±u ) 

(6) 

where Ad is the anisotropic term of the hyperfine cou- 
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FIG. 7: Temperature dependence of the ratio D/loh between 
the spin diffusion coefficient D and the exchange frequency 
ljh in [YPC2] compound as derived from the slopes in the 
1/Ti vs (uj /2-n-y 1/2 plots in Fig. \E\ The solid line gives the 
best fit according to D oc exp(A/T) with A = 1.2 ± 0.4 K. 

pling, which hereafter shall be neglected since A 2 ^> A\. 
Then just the second term in square bracket can be con- 
sidered. In Eq|5] xo is the static uniform susceptibility 
per spin and oj e = wo7e/7 is the electron resonance fre- 
quency. This means that during the nuclear relaxation 
process a simultaneous flip of the electron and nuclear 
spins occur, involving an energy exchange h(uj e ± (jJq), 
and 1 /Ti thus probes the spin excitations at a frequency 
close to oj e . 

In a one dimensional system, the spectral density at 
uj e is characterized by a low-frequency divergence given 
by!2 



, 1 ,U) C + y/uf+UJ* 1/2 
Je > - 9~, 5 ) 



(7) 



where w c is a low-frequency cutoff accounting for the fi- 
nite spin anisotropy and/or inter-chain coupling, while 
D is the spin diffusion rate. In Fig. (|6]), the 1/Ti 

— 1/2 

is plotted as a function of v . The observed linear 
trend further proves the one-dimensional nature of the 
antiferromagnetic correlations. Moreover, the absence 
of a low-frequency flattening in 1/Ti plot indicates that 
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spin diffusion occurs in the electronic frequency range 
uj c <C uj e <C D. Thus, from the slopes of the curves it is 
possible to deduce the spin diffusion coefficient at differ- 
ent temperatures (Fig. [7]) considering A ~ 180 Gauss and 
neglecting uj c <C w e in Eqs <[6]{7l) . The estimated spin dif- 
fusion coefficient is of the order of the exchange frequency 
loh and it is found to progressively decrease with tem- 
perature and to become nearly constant above 20 K. It 
is interesting to observe that D oc exp(A/T), namely the 
behaviour expected for onc-dimcnsional antifcrromagnets 
in the presence of a spin-gap A between singlet and 
triplet excitations. 20 Here we find that A = 1.2 ± 0.4 
K suggesting that a small gap, either due to compet- 
ing exchange interactions or to a dimerization might be 
present in [YPC2] . It is interesting to observe that, at 
low temperature, when the Zeeman energy Huj e ~ A the 
breakdown of Eq[7] is noticed. In fact, in Fig. ([BJ one 
clearly notices that at T — 1.6 K the linear behaviour 
is no longer obeyed at high fields (i.e. low values for 
\/2it /ujq) and 1/Ti ceases to decrease with increasing 
field. This could be due to the modifications in the spin 
correlations induced by the magnetic field for fvjj e ~ A, 
possibly associated with the progressive closure of the 



spin gap. 

In conclusion, from magnetization, X H NMR param- 
agnetic shift and T\ measurements we have derived the 
magnitude of the antiferromagnetic exchange interaction 
in [YPC2] compound and found an overall good agree- 
ment. The low-energy spin excitations are of diffusive 
character and characteristic of one-dimensional antiferro- 
magnets. From the temperature dependence of the spin 
diffusion rate derived from 1/Ti vs. H measurements it 
was found that a spin-gap around 1 K might be present 
in this compound. 
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